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1. Introduction 

The masses of black holes in the nuclei of galaxies are measured by observing how they accel- 
erate test masses in their vicinity, as is always the case in astronomy. The "test masses" in the case 
of galactic nuclei are stars or gas in the nuclear regions. The particular advantage of using stars 
as the test masses around black holes is that they are subject only to gravitational forces. Unfortu- 
nately, stellar dynamics can be modeled accurately only with high angular-resolution spectroscopy; 
the black hole radius of influence, 

R BU = 2GM B hM 2 , (1.1) 

where a* is the stellar velocity dispersion of the host galaxy bulge, must be resolved, or at least 
nearly so. Measurements of gas motions, on the other hand, allow us to probe much closer to 
the black hole, on scales far below the black hole radius of influence in the case of reverberation 



mapping [ ]15| , [76Q . However, gas also responds to non-gravitational forces, such as hydromagnetic 
acceleration and radiation pressure, which must be accounted for. 

Methods that rely on measurements of the motions of stars or gas that are accelerated by 
the black hole are "direct methods." These include modeling of both stellar and gas dynamics 
and reverberation mapping. We also employ "indirect methods" of mass determination that mea- 
sure observables that are correlated with black hole mass: these include the relationships between 



black hole mass and stellar bulge velocity dispersion (the Mbh~<7* relationship [ p4 , 38L p4|]), black 
hole mass and host-galaxy bulge luminosity (the MBH-^buige relationship, which is also sometimes 
known as the Magorrian [ ]6T| ] relationship), the fundamental plane, and AGN-specific relation- 
ships, such as that between the AGN luminosity and the size of the broad-line region (BLR), as 
discussed below and by Bentz []|]. It is also useful to distinguish among "primary," "secondary," 
and "tertiary" methods of mass determination: primary methods are those that require the fewest 
assumptions. Certainly, in the case of supermassive black holes, mass measurements based on the 
proper motions and radial velocities (Sgr A*) or megamasers (NGC 4258) are primary. Also gen- 
erally regarded as primary in the measurement of nuclear supermassive black holes are stellar and 
gas dynamics. Reverberation mapping, however, as it is currently practiced, requires an external 
calibration of the zero-point for its mass scale through correlations between the black hole mass and 
host-galaxy properties that are assumed to be the same for the host galaxies of active and quiescent 
black holes; for this reason reverberation mapping is a "secondary method," although it is still a 
"direct method." However, as discussed below, reverberation mapping has the potential to become 
a primary method. Since it will be the only method that is potentially extendable to high redshift 
and the method most applicable to NLS Is, it will be the main focus of the rest of this contribution. 

2. Results from Reverberation Mapping 

The desired outcome of reverberation mapping is usually expressed in terms of a "velocity- 



delay map" (also called the "transfer function" [15]) *^(t, Vlos) which is the six-dimensional phase 



space of the BLR (velocity field and geometry) projected into the two observables, time delay % 
and line-of-sight (i.e., Doppler) velocity Vlos- The data requirements for successful recovery of 
*P(t,Vlos) from spectrophotometric monitoring are fairly stringent [[f7|] and it is therefore only 
recently that high-quality maps are starting to appear fll4|, 53l. It is much more common to extract 
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only the mean response time for an emission line integrated over all Doppler velocities. Simple 
cross-correlation of the continuum and emission-line light curves gives a mean response time for 
a particular emission line, and this mean response time is weighted toward physical conditions 
where both the emissivity and responsivity (marginal change in emissivity in response to a contin- 
uum change) are high. Even in this simplified analysis, advances continue to be made as improved 
mathematical descriptions of AGN variability [53, ^(J have led to more observationally motivated 
statistical modeling of the continuum behavior in the gaps in coverage in the real time-series sam- 



pling. Recent improvements in this methodology [104] allow self-consistent determination of the 
uncertainties and, since multiple lines with different lags can be dealt with simultaneously, make 
it possible to back-fill some of the gaps in time series. The new methodology yields lags that are 
largely consistent with those based on current cross-correlation practices, but with reduced uncer- 
tainties in the lag measurements and, indeed, more robust estimates of these uncertainties. 

At the time of this workshop, reverberation lag measurements have been made for around 50 
AGNs. Since the compilation of Peterson et al. [|82|], there have been new reverberation results 
from MDM Observatory and Crimean Astrophysical Observatory [0, ||, ||, |7[ ||, H El> ES' 
the Lick AGN Monitoring Program (LAMP) [§, jjj, [l3|, 0, |] and other programs ||| |9l]] that 



have concentrated primarily on the B aimer lines in low-redshift AGNs. There have also been new 



measurements of the C IV A 1549 response in the dwarf Seyfert NGC 4395 [83] and in one high 



redshift quasar [52]. The Mg II A 2798 response has been measured reliably in only one case, NGC 



4151 [69] 



The major conclusions that can be drawn from these studies are: 

1. Different emission lines respond on different time scales, with lines characteristic of the 
highest ionization gas responding first, providing evidence for ionization stratification of the 
BLR. The size of the BLR, as measured by a particular emission-line time delay T, is inferred 
to be Rblr = ex. 



2. 



In every case where the time delays for multiple broad emission lines have been measured 
[77, 78], 58, 13] there is a tight anticorrelation between line width AV and time lag T such that 
the product AV 2 T is constant (Figure 1), as expected if the dynamics of the BLR is dominated 
by the gravitational potential of the central black hole. Strictly speaking, it means that the 
BLR gas is responding to an inverse-square force or forces, so radiation pressure might 
also play a role [ ]62| ] (see Marconi's contribution to these proceedings for a more thorough 
discussion [|63|]). 



3. The time lag for the emission-line response to continuum variations is longer in more lumi- 



nous objects Q50t , plQ . After separating the AGN luminosity Lagn from the starlight contam- 
ination by the host galaxy [^, |TI|], the relationship between these quantities is consistent with 

1 /2 

^blr x ^agn- This * s well-established only for Hj8, but the more limited reverberation results 



for C IV [82, 831, p2J are consistent with the same relationship. This radius-luminosity (R-L) 



relationship provides the underpinning for indirect methods of measuring AGN black-hole 
masses over cosmic time [^|, 95]. 



4. In every case in which the response of a particular emission line has been measured on 



multiple occasions, the timescale for response of the emission line is consistent with T°=L 



1/2 
UV 
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Figure 1: Line-width vs. emission-line lag (in days) for four reverberation-mapped AGNs. In each case, the 
data are consistent with a virial relationship, i.e., V T~ 1//2 . The dotted lines are show the best-fit slopes 
and the solid lines show the best-fit virial relationship [fnl [78|. 



where Luv is the luminosity in the observable UV, as close to the ionizing continuum as 
possible [Kl[]. The emission-line width also changes in a way that is at least approximately 
consistent with a constant value of the product AV 2 T. 

It is worth mentioning that the strong optical Fe II blends that are prominent features of Seyfert 1 
spectra are seen to vary over time scales that are long compared to reverberation time scales, but 
do not appear to vary much if at all on reverberation time scales 57]. This may be because the 



emissivity and responsivity distributions for these lines are not highly localized, as they seem to be 
in the case of most of the prominent emission lines in AGNs. In other words, these features vary, 
but they do not reverberate. 

Results (2) and (4) above suggest that we can measure the central black hole mass by combin- 
ing the emission-line width AV and time lag T, i.e., 

'AV 2 ct s 



M 



BH 



f 



G 



(2-1) 
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The quantity in parentheses is based on the observable quantities and has units of mass; we refer 
to it as the "virial product." The factor / is a dimensionless number of order unity into which 
we subsume all of our ignorance of about the geometry and kinematics of the BLR; what we are 
attempting to do is to characterize a probably complex BLR with two observables that we hope 
will somehow average over the complexity well enough to yield a reasonable mass estimate. Let's 
consider the elements of this equation individually: 

• As noted earlier, T is the mean time lag for a given emission line. Practically speaking, 
it is usually the centroid of the continuum/emission-line cross-correlation function. It is 
unbiased as a function of BLR inclination as long as the BLR is axisymmetric and emits 
isotropically (and even when it does not emit isotropically, the correction factor is relatively 
small compared to other projection effects). 

• The line width AV is more problematic, as there are multiple issues to consider. The first 
of these is what line-width measure to use? Commonly used measures are (a) full-width at 
half maximum (FWHM) and (b) line dispersion G\i ne , the second moment of the line profile. 
Both measures have liabilities: FWHM is more sensitive to both random noise, unremoved 
narrow-line components, and absorption due to intervening gas (at least in the case of reso- 
nance lines like C IV A 1549) and aii ne is more sensitive to blending with other features 1 . Sec- 
ond, given an array of spectra from a reverberation-mapping monitoring campaign, we have 
two obvious alternatives for forming the spectrum from which we measure the line widths. 
The first of these is the mean spectrum, which we construct by determining the average flux 
in each wavelength bin of the spectrum over the duration of the monitoring campaign. The 
second of these is an "rms spectrum" that is comprised of the root-mean-square fluxes in 
each wavelength bin of the spectrum over the duration of the campaign. The rms spectrum 
affords two particular (related) advantages: (a) the rms spectrum isolates the emission-line 
gas that is actually responding to the continuum variations and omits the unchanging parts 
(thus automatically excluding narrow-line components and host-galaxy starlight), and (b) we 
find in practice that the variable parts of the emission lines are much less blended with other 
features, mostly because some features such as the optical Fe II lines do not vary on reverber- 
ation time scales, but partly because the widths of the emission lines are typically narrower 
in the rms spectrum than in the mean spectrum (possibly because the very highest velocity 
BLR gas is optically thin [[33|, |^]). A nice comparison of a mean and an rms spectrum is 



shown in Grier et al.'s contribution to these proceedings [}43[]. 



• The scaling factor / converts the convenient virial product, based on the observables, to an 
actual mass. If the BLR is a flat Keplerian disk observed at inclination i, the scaling factor 
would include a 1/ sin 2 i term to account for the velocity projection. It could also include a 
correction factor for z if the line emission is anisotropic. Also — and this is a critical point 
that is often misunderstood — the value of / depends on what line-width measure is being 
used. This is an extremely important point that I will return to below. 



We are currently experimenting with using interpercentile widths [IOC] as an alternative characterization of line 



width, as have Fine et al. [poj], but it is premature to conclude anything at this point. 
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Table 1: Black Hole Masses (Units of 10 6 M e ) 



Galaxy 


NGC 3227 


NGC 4151 


Direct methods: 
Stellar dynamics 
Gas dynamics 
Reverberation 


7-20 [23 
20+}° [At 
7.63 ±1.7 [ 


] 

] 

29 


<70 [73] 
30+7^ [(§] 
] 46±5 0] 


Indirect methods: 

Mbh-ct* 
R-L scaling 


25 
15 




6.1 
65 



It should be clear from this discussion that the scale factor / is different for every AGN; indeed, 
in principle it could be somewhat different for different emission lines in the same AGN, but not 
too much different since, as shown in Figure 1, the virial product is approximately constant for all 
emission lines. Without additional information on the kinematics, geometry, and inclination of the 
AGN, we cannot determine / for a particular source. 

We could also determine / if we had an independent measure of the black hole mass. Un- 
fortunately, there are very few cases where alternative direct measurements exist (see below), and 
none of these are high precision measurements. We can, however, use other indirect methods to 
estimate the black hole masses in reverberation-mapped AGNs and use these to determine the scale 
factor. However, since the fidelity of individual mass measurements by indirect methods is low, 
we can realistically hope to obtain only a mean value (/) for the reverberation-mapped sample and 
thus effect a statistical scale factor for the reverberation-mapped sample. This has been done by 
assuming that the Mbh~<7* relationship for AGNs is the same as it is in quiescent galaxies. It is 
well-known that reverberation-based black hole masses scale with stellar bulge velocity dispersion 
in a fashion similar to that seen in quiescent galaxies [39, 35, 72, 70, |103 ] so it is not too much of 
a stretch to assume that the zero-points for the AGN and quiescent-galaxy are identical. The most 
recent analysis gives a value of (/) = 5.25 ± 1.21 [103] (though see [ |40| ] and §43). The scatter 
around the AGN Mbh~<7* relationship is estimated to be ~ 0.4dex, which is thus an estimate of 
the typical uncertainties in reverberation-based masses, assuming that the intrinsic scatter in the 
Mbh~<7* relationship is much smaller than this (but more on this later). 



The reverberation-based black hole masses in the literature are computed from equation (2.1) 
with AV taken to be the line dispersion in the rms spectrum, X is the cross-correlation centroid, and 
/ is the value based on the Mbh~<7* relationship [72, 103]. In a small number of cases, we can com- 
pare the reverberation-based masses with those from other direct methods. Two cases are shown in 
Table 1, NGC 3227 and NGC 4151. In both cases, the direct methods are in reasonable agreement, 
particularly considering the factor of three or so systematic uncertainty in the reverberation-based 
masses and probably comparable systematic uncertainties in the other primary methods. Another 
check on the accuracy of reverberation-based masses is to compare the MBH-^buige relationship for 
AGNs to that for quiescent galaxies [10]. The agreement is found to be very good. 

Use of a mean value (/) averages over the specific values of / for individual sources: if it is 
unbiased, (/) will result in equal numbers of mass overestimates and underestimates, but these will 
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not be a function of either time delay or line width. 



3. Indirect Mass Measurements Anchored by Reverberation Mapping 

As noted above (result 3 in the last section), reverberation mapping has revealed a tight rela- 
tionship between the AGN continuum luminosity and the radius of the BLR ^ of the form 

^BLR^AGN; (3-1) 

where empirically we find a 0.5, which is, as is often noted, consistent with the theoretical 



expectation. In the early days of photoionization equilibrium modeling of the BLR, Davidson [21] 
and others began to recognize that photoionization models could be parameterized by (a) the shape 
of the ionizing continuum, (b) the elemental abundances, (c) the particle density «h and column 
density of the photoionized gas, and (d) an ionization parameter defined by 

U= m (32) 

where R is the separation between the source of ionizing photons and the photoionized gas (i.e., the 
BLR radius) and 2(H) is the number of ionizing photons produced per second by the photoionizing 
source, i.e., 



2(H) = f ^dv (3.3) 



and the integral is over all hydrogen-ionizing frequencies. To some low-order of approximation, 
AGN spectra are self-similar over many orders of magnitude in luminosity, which suggests that the 
BLRs of all AGNs are characterized by similar values of U and «h- This leads to the prediction 
that 

( 6(H) \ VZ x L i/2 

if we can use Lagn> the luminosity at some observable wavelength, as a proxy for the ionizing 
luminosity. We characterize this expectation as "naive" since some of the implicit assumptions 
are quite approximate: we know that the shape of the ionizing continuum should change with 
black hole mass and that the BLR is complex and cannot be characterized by fixed values of U 



1/2 



and «h M, 55]. In the early photoionization equilibrium calculations in the 1970s, the size of 



the BLR was not considered to be an important parameter as the models suggested that the BLR 



sizes in quasars would range from one to a thousand light years or so [ |56| , |22| ], which would be 
unresolvable in distant quasars in any event. I suspect that Davidson does not get enough credit 
for his insight on the ionization parameter and the R-L scaling relationship that follows from it 
because he wrote the ionization parameter in terms of the ionizing flux (U oc Fi on /nu) rather than 
the ionizing luminosity (U °< Li on /R 2 nn), as it was written later [67], so the BLR radius did not 



appear explicitly. Nevertheless, the prediction of the R-L relationship was so well-known that 



a decade later it was often quoted without attribution (e.g., [64]) when the BLR size became an 



interesting parameter because the first emission-line variability studies suggested that the BLR was 



much smaller than predicted by the photoionization equilibrium models [79]. The R-L relationship 



was looked for as a prediction of photoionization theory even when the first crude reverberation 
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lags were measured [ 54 ] ; it was always understood, of course, that the central mass could be easily 



estimated from a single quasar spectrum via equation ( |2~ll ) by using /?blr from the R-L relationship 



and AV from an emission-line width measurement. This "photoionization method" J98| , pQp thus 
provides a significant shortcut to estimating AGN black hole masses 2 . The general methodology 
has been extended to other emission lines, in particular, C IV X 1549 |94], and Mg II A2798 
[]68|]. Similar R-L relationships using various proxies for the continuum luminosity have also been 
shown to be viable J4l|]. 

The bottom rows of Table 1 show two indirect measurements of the black hole masses in NGC 
3227 and NGC 4151. The indirect methods are in good agreement with the direct methods in the 
case of NGC 3227, but not quite as good in the case of NGC 4151. But generally the picture holds 
together to an accuracy of ~ 0.5 dex. 

4. Implications for NLSls 

4.1 Inclination or Eddington Ratio? 

There are two competing explanations for the properties of NLSls: 

1. NLSls might be low-inclination sources, viewed nearly pole-on. If the broad-line emission 
arises primarily in a rotating Keplerian disk, for example, the projected rotational speed is 
decreased by a factor of sin i, thus accounting for the unusually narrow broad lines. 

2. NLSls might be AGNs that are accreting at relatively high Eddington ratio: for a given 
luminosity, objects with the narrowest emission lines have the lowest black hole mass. 

There is certainly evidence that in general broad lines widths are affected by inclination. The 
widths of the B aimer lines in core-dominated radio sources (those at low inclination) are systemati- 
cally lower than those in the spectra of lobe-dominated (high-inclination) sources [ ]102j ]. Similarly, 
radio sources with flat spectral indices (low inclination sources) have narrower line widths than 



radio sources with steep spectral indices Q49p . Of course, the fact that "broad" lines are never nar- 
rower than ~ 1000km s~ ! means that there must an an axial component to the BLR velocity field, 



but it still seems that the dominant motion might be rotational. Decarli et al. [ |24| ] make a good case 
that such a model can explain many of the properties of NLS Is, although I don't believe inclination 
can explain everything. Also, absorption-line and reverberation studies [19, both suggest that 



the well-known NLS1 NGC 4051 is indeed observed at low-inclination. And, very importantly, 
some NLSls have now been detected in very high energies by Fermi (QlL |37[ ], suggesting that we 
are looking down the axis of these systems. 

Is inclination sufficient to account for the NLS 1 phenomenon? Probably not, in my opinion. 
I'll cite two specific examples. One of the few radio-loud reverberation mapped AGNs is 3C 120, 



2 This method for estimating quasar masses has been sometimes, misleadingly in my opinion, referred to as the 
"Dibai method" Jig], Dibai noted that since AGN emission lines have the same equivalent widths regardless of 
luminosity L\ i ne °z ^continum • one then assumes constant emission-line emissivity per unit volume £, then L\[^ — 
£47r7?g LR /3, or Rblr 01 ^ccmt anc ' ^bh x ^^ 2 ^cont- So while the predicted sizes for the BLRs in nearby AGNs were 
reasonable, the physics is incorrect and the functional form of the R-L relationship is wrong. Certainly Dibai did a lot 
of important work that was admittedly often overlooked in the West, but this incorrect attribution overlooks the more 
correct physical insights provided by Kris Davidson, Gordon MacAlpine, Chris McKee, Bruce Tarter, and others. 
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Table 2: Black Hole Mass Measurements for Mrk 1 10 



Method 


Mass (Units of 1O 6 M ) 


Reverberation: 


25±6 


Mbh-ct*: 


4.8 


Gravitational redshift: 


14±3 



which has a relativistic jet with an inclination i < 20°. The inclination correction to the mass 
is at least a factor of 10 for this inclination; but 3C 120 nevertheless falls right on the Mbh-<7* 
relationship rather than below it. Mrk 1 10 is another interesting case, an NLS1 with an independent 
black hole mass measurement based on the gravitational redshift of the emission lines [|^]. As seen 
in Table 2, the reverberation mass is the largest of the three available mass estimates: we would 
expect it to be the smallest if the narrowness of the B aimer lines is attributable to inclination only 
as the other two measures should be independent of inclination. 

Consider now the alternative explanation, that the NLS1 phenomenon is a manifestation of 
low black-hole mass at a given luminosity. Figure 2a shows the mass-luminosity relationship for 
reverberation-mapped AGNs, with the NLSls highlighted. The NLSls are generally high Edding- 
ton ratio (i.e, accretion rate relative to the Eddington value) objects, but they are also generally 
low-luminosity sources. This, we shall see, is an accident of how the NLS 1 class was defined. 

Let's suppose that the common physics of the NLS 1 phenomenon is captured by the Eddington 

1 /2 

ratio. The H/3 line width is given by the virial equation AV °< (Mbh/^blr) ■ From the R-L 
relationship of the last section and the definition of the Eddington ratio m = M/M^dd x M/Mrh, 
we see that the dependence of the line width is 



... lM m \" 2 /MbhV' 2 /MbhV' 4 



"AGN, 



What this tells us is that for a fixed Eddington ratio m, the emission lines are broader for more 
massive black holes. If this is correct, then we should be modifying our definition somewhat, 
since higher mass black holes will have broader lines even though the Eddington ratio is fixed 



4, 32]. Rather than modify the NLS 1 definition, I'll just define (arbitrarily) a larger class of "high 



Eddington ratio sources (HERS)" as those with 

Ay ^(T^) 1/4 2000kms "- (4 " 2) 

The HERS are highlighted in Figure 2b. If the NLSls are a high Eddington ratio phenomenon, 
then they are the low-luminosity subset of the larger HERS class. But taken as a group, they do not 
stand out from the other objects in any properties except line width and profile and the parameters 
derived directly from these parameters; this is consistent with Bentz but rather at odds with the 
results of other investigations (see §|4~3|). 



The bottom line is probably that NLSls are yet another "mixed bag" (cf. [101]) of sources 
that include both high Eddington rate accretors and low-inclination AGNs. Inclination effects are 
clearly present, as are some NLS1 characteristics that cannot be explained away as inclination 
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Bolometric luminosity 




Optical luminosity XL X (5100A) 

Figure 2: Mbh-^agn for the reverberation-mapped AGNs. Lines of constant Eddington ratio are shown 
for m = 0.01, m = 0.10 and m = 1.0. Multiple Lagn values for a fixed Mbh indicate multiple reverberation 
results for some sources (notably NGC 5548 at ~ 7 x 10 7 M fl ). The upper scale shows bolometric luminosity 
based on a nominal correction Lbol = 9ALagn- In the upper panel, NLSls are shown as open cirles. In the 
lower panel, high Eddington rate sources (HERS) are shown as open circles. 



effects. Indeed, the case can be made that many NLSls are actually both pole-on and have high 



Eddington ratios [17] 



4.2 Characterizing the Line Width 

So far we have side-stepped the question of what is the better measure of the line width, 
FWHM or aii n e? In this case, "better" means which measure, if either, gives a black hole mass that 
is unbiased. It is inescapable that at least one of these measures will introduce a bias into the mass 
scale because the shape of AGN emission lines is a function of their width. This is illustrated in 
Figure 3, which shows FWHM/aii ne as a function of FWHM, where again we have highlighted the 
HERS. Clearly the choice of which measure we use for computing the masses is important. This 
is illustrated very clearly in Figure 4, which shows a direct comparison of black hole masses based 
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Figure 3: The line-width ratio FWHM/ah ne as a function of FWHM for the Hj3 line in the rms (top) and 
mean (bottom) spectra of reverberation-mapped AGNs. In both panels, a dotted line shows the value for a 
Gaussian profile, FWHM/ Oji ne = 2.35. HERS are shown as open circles. 



on FWHM with those based on aii ne ; the HERS have smaller masses for FWHM than for a une and 
vice versa for the remainder of the AGNs. 

An argument can be made that a une is the less-biased line-width measure, at least in the case 
of lines measured in rms spectra [82, 18], the variable part of the spectrum. However, in individual 
spectra of AGNs, generally the line wings are badly blended with other features, and introduces 
systematic errors that are a function of line width [j26|], but the errors introduced are still smaller 
than the systematic differences that arise from using FWHM instead of aii ne - 

The clear trend illustrated in Figure 3 demonstrates that it is possible to correct FWHM to 
aii n e ; this was, in fact, explicitly proposed by Collin et al. [18]. There are other formulations that 
have suggested that the dependence on line width in equation (2A) should be written as FWHM 7 , 
where 7 is a free parameter, which in practice turns out to have a value slightly smaller than two 
]; in essence, this is making the same correction between FWHM and a une , since the latter is 



what was used to compute the reverberation-based masses [82]. 

It is important to sort this out as incorrect parameterizations will give misleading results. While 
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M Ba (o lhl -based) 

Figure 4: A direct comparison of black hole masses based on different line-width measures, FWHM and 
the line dispersion <7ij ne , in both cases based on measurement of the Hj3 line in the rms spectrum. HERS 
are shown as open circles. Relative to Cii ne , FWHM yields lower masses for HERS because of the higher 
kurtosis of the Balmer-line profiles in HERS spectra. 



the issue has not been completely resolved at this point, there is a current controversy that in any 
case illustrates the point. Steinhardt & Elvis [92] show that if one plots black-hole mass versus 
AGN luminosity for the SDSS DR5 quasars, in each redshift bin it is found that while the lower- 
mass AGNs approach the Eddington limit, the higher mass AGNs do not. In other words, the 
axis of the distribution of quasars on the Lagn-^bh diagram (similar to Figure 2, with the axes 
interchanged), does not run parallel to lines of constant m, but is rather more shallow. The absence 
of higher mass AGNs with high values of the Eddington ratio is referred to as the "sub-Eddington 
boundary." The mass measurements used by Steinhardt & Elvis are from Shen et al. [9C], in 
which the emission lines were fit with a pair of Gaussians, one for the narrow component and one 
for the broad component, and the black hole mass was computed using equation (2T) with the 
FWHM measurement for the broad-component best-fit Gaussian. In an independent study, Rafiee 
& Hall [ |85l ] went back to the original DR5 spectra and measured Oii ne from the original data. 
When they use aii ne to compute the black hole masses, they find that the sub-Eddington boundary 
disappears, and that in each redshift range the distribution of AGNs in the Lagn-^bh diagram is 
along constant values of rh. This, of course, does not prove that the sub-Eddington boundary is an 
artifact of how the line widths were characterized, but it is certainly suggestive since the Rafiee & 
Hall result seems intuitively more likely to be correct. Indeed, Figures 3 and 4 show that using 
FWHM instead of aii ne will certainly lead to overestimating the masses of AGNs with the broadest 
lines and underestimating the masses of the AGNs with the narrowest lines, which would produce 
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exactly the sub-Eddington boundary phenomenon. A Gaussian fit (rarely a good description of a 
broad-line profile) will have the same effect (see Figure 3), so it's not immediately clear whether 
the problem is principally the Gaussian assumption or the use of FWHM. 
Is the sub-Eddington boundary real? Perhaps. But I doubt it. 



4.3 Host-Galaxy Properties 

Currently the reverberation-based mass scale relies heavily on correlations between black hole 
mass and host-galaxy properties. It is explicitly assumed that these relationships are the same in 
active and quiescent galaxies. A particular concern is the nature of the Mbh~<7* relationship, which 
currently provides the zero-point calibration for the reverberation mass scale. Batcheldor [Q] has 
argued that the Mbh-<7* relationship is a selection effect caused by preferentially observing sources 
with a resolvable radius of influence (equation though Gultekin et al. p5| ] disagree. Indeed, 
the existence of an AGN Mbh~<7* relationship also argues against this as a selection effect, since 
there is only a subtle radius-of-influence bias in the AGN sample 3 . 

There is also active discussion about the host galaxies of NLS Is and whether or not there are 
differences between these and the hosts of other AGNs. Crenshaw, Kraemer, & Gabel | ^0| ] claimed 
that NLS1 hosts are disproportionately barred galaxies, and larger samples seem to confirm this 
J7T1]. This raises concerns as there is evidence [[[(]] that the Mbh-c* relationship is different for 
barred galaxies than for non-barred and elliptical galaxies. There are also indications that the hosts 
of NLS 1 s have pseudo-bulges rather than classical bulges [ 66 , 74 , |65 , 75 ] , and pseudo-bulges might 



[]5|] or might not [ |48| ] not follow the Mbh~<7* relationship. This is a potential problem not only for 
the overall reverberation-mapping mass scale, but for the relative masses of NLS Is compared to 
other AGNs. 



4.4 Toward Higher Precision Masses 

At the present time, the masses that can be obtained from single-epoch spectra are probably 
accurate to about the 0.5 dex level, but we start to see hints of bias in the mass scale at this point. 
Proper characterization of line widths in individual (single-epoch) spectra is probably the single 
most significant problem, although this does not seem to be generally appreciated. Beyond this, 
there are several important questions that need to be addressed as we strive to make better black hole 
mass measurements, and some of these are addressed in other contributions to these proceedings: 

• Are black hole mass/host-galaxy relationships the same in AGNs and quiescent galaxies? Do 
AGNs and quiescent galaxies differ in incidence of pseudo-bulges, and how does this relate 
to black hole growth and star-formation rates [ [86] , |87|]? 

• Are black hole mass/host-galaxy relationships the same in NLS Is (or HERS) as in other 
AGNs? 

• How reliable are the scaling relationships over luminosity and redshift? 

3 Specifically, the most luminous objects in the reverberation database were selected at least in part for their apparent 
brightness, so there is a bias toward HERS, as is well-known for the Palomar-Green sample from which these are drawn. 
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• Does failure to account for radiation pressure lead us to underestimate Mbh? NLSls (or 
HERS) are the likely testing ground for this because the radiation pressure term in the mass 
equation J63| ] will be relatively highest in these sources. 



What problems do we encounter by using different emission lines to determine Mbh? In 
particular, there are lingering doubts about the use of C IV, despite good arguments that it 



is generally a good mass-indicator [30, My own observation (see also [30, ^]) is that 
much of the evidence that C IV is problematic is based on data that are too poor to address 
the issue. 

• How can we resolve the inclination/Eddington ratio ambiguity in NLSls (or HERS)? 
Which brings us . . . 

5. Back to Reverberation Mapping 

While we expect that ongoing research will lead toward higher accuracy and higher precision 
black hole masses, improvements in reverberation mapping is the direct path to address all of our 



concerns about black hole masses. The velocity-delay map for the Balmer lines in Arp 151 [14] 



indicates a rich, complex structure and a preliminary velocity-delay map for the NLS 1 NGC 405 1 



[30] suggests that the Balmer-line emitting region in this source is a nearly face-on disk. There 



are now several data sets of comparable quality []43[], so additional velocity-delay maps can be 
expected in the near future. As more and better velocity-delay maps become available, it will be 
possible to model the kinematics and geometry of the BLR directly from the reverberation data and 
thus determine the central black hole masses without recourse to other methods; this is what will 
be required for reverberation mapping to become a primary method. 

The problem with reverberation mapping is, of course, simply that it is both expensive (in 
terms of telescope time) and risky (since the pattern of variability in a given source is unpre- 
dictable). In the future, it will probably be possible to employ more clever strategies for rever- 
beration campaigns, for example, using photometric monitoring from synoptic surveys like LSST 
to trigger campaigns at opportune times. And given good photometric coverage, spectroscopic 
monitoring can be sparser and tailored to the goals of the program (for a good example of an 
opportunistic program, see 

But it remains true that reverberation mapping is currently the only way to measure the masses 
of black holes in galaxies out to cosmologically interesing distances and the only way to probe the 
low-mass end of the distribution, as illustrated in Figure 5. Even with a diffraction-limited 30-m 
telescope at //-band, it is impossible to resolve the black hole radius of influence for Mbh < 10 8 M Q 
at the distance of the Coma cluster. Reverberation mapping, by substituting time resolution for 
angular resolution, can extend direct measurements of black hole masses to much larger distances. 
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D A (Mpc) 

Figure 5: Mbh versus luminosity distance for nuclear black holes with direct measurements of their masses. 
Reverberation-based masses are shown as open circles and masses measured by stellar or gas dynamics are 
shown as filled circles. The limiting black hole radius of influence as a function of mass is shown for the 
cases of Hubble Space Telescope at the Call triplet (~ 8600 A), an 8-m telescope at //-band, and a 30-m 
telescope at //-band. A bare handful of AGNs are close enough to resolve the radius of influence and thus 
measure their black hole masses by the same techniques used on relatively local galaxies. 
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